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^ , Abstract: Qualitative features of the mean-field theory of superconductivity in a strongly 

disordered system of fermions with short-range attraction are discussed. In this limit the 
effective theory at low energies is entirely bosonic, and I consider both the artificial infinite- 
range limit and the more realistic case of "nearest-neighbor" hopping of bosons between the 
'^ ■ localized states. In the infinite-range case the mean-field theory is exact, and the supercon- 

Q . ducting gap is uniform in space. There is a smooth BCS-BEC crossover with decrease in 

CJ I density, at weak enough disorder; at moderate densities, or larger disorder, the mean-field 

ground state is the BCS-like locahzed superconductor. In the latter case, the gap is highly 
^ ■ non- uniform in space, but surprisingly is finite everywhere at T < Tmf- I find that the mean- 

ly-^ ■ field transition temperature Tmf > always, and argue that the superconductor-insulator 

O , transition at T = in models with net attraction between fermions is in the universality 

^ ! class of "dirty-bosons". 
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1 Introduction 

The effect of disorder on superconductivity has been one of the fundamental issues in the 
field, which despite a lot of effort has still not been completely understood. It is well known 
that weak non-magnetic impurities do not affect much the transition temperature of a s- 
wave superconductor |jl|, p|, since one only needs to replace the plane- waves in the BCS 
variational ground state with the pairs of time-reversed exact eigenstates of the random 
single-particle part of Hamiltonian, to essentially preserve the result of the BCS theory. The 
eventual breakdown of the Anderson theorem with increase of disorder was studied in mid- 
eighties by several groups 0, Q], [|]. It was found that the BCS variational ground state 
continues to have lower energy than the non-superconducting state until the extreme limit of 
site-localization is reached, when the system becomes a gapless insulator of localized, tightly 
bound, singlet electron pairs. The inclusion of quantum fluctuations should lead to loss of 
long-range order in the ground state at a finite value of disorder, and the low-energy sector 
of the theory was argued to be equivalent to the system of disordered bosons with hard-core 
repulsion [Q . The question of critical behavior at the quantum (T = 0) superfiuid-Bose glass 
transition in the system of dirty bosons has been of great interest [^], since it also represents a 
paradigmatic case of poorly understood interplay of localization and interactions. Recently, 
a systematic theory of the superfiuid-Bose glass quantum critical point was developed based 
on the expansion around the lower critical dimension d = 1 for the transition ^ . 

Despite the success p of the bosonic theory of the quantum superconductor-insulator 
transition in naturally accounting for the observed scaling of resistivity and for the apparent 
universality of the critical conductivity, there is no consensus yet on the basic issue of its 
physical correctness. The bosonic theory assumes that electron pairs do not dissociate at the 
point of transition, but that it is only the long-range order in the phase of the superconducting 
order parameter that disappears. This is exactly what happens, in the strict renormalization 
group sense, at a finite temperature superconducting transition. A complementary view 0, 
supported by a different set of experiments [l^, is that disorder enhances the effect of 



Coulomb interaction between electrons, to the point at which the net interaction becomes 
repulsive and the pairs break apart. Which of the two pictures should apply to a given 
experimental situation is at the present time far from being clear. 

The purpose of the present paper is to clarify some of the conceptual issues behind 
the dirty-boson theory of the superconductor-insulator transition. In particular, I study the 
mean-field transition temperature in a model of strongly disordered fermions with a fixed net 
short-range attraction. Such a simplified model should presumably always show a quantum 
phase transition from a superfiuid into a Bose glass of electron pairs, at T = and with in- 
crease of disorder. One expects that the mean-field transition temperature, which should be 
understood as a crossover scale where the amplitude of the superfiuid order-parameter starts 
forming, would remain finite at the critical value of disorder, i. e. that the superconducting 
gap does not collapse at the same point where the phase-coherence disappears (see Figure 1). 
This point is studied here by revisiting the mean-field theory of a disordered superconductor 
M, but this time considering the problem at finite temperatures. Can the superconducting 



mean-field transition temperature Tmf in a disordered Fermi system with purely attractive 
interaction between fermions vanish? Naively, one would allow this possibility at low densi- 
ties of electrons, where the hard-core repulsion between bosonic electron pairs becomes less 
important, and the system should behave similarly to the gas of nearly independent bosons 
moving in a random external potential. The Tmf would then correspond to the temperature 
where the Bose-Einstein condensation (BEC) into the lowest single-particle state takes place 
||TT| . In this dilute limit, if the density of random eigenstates at the chemical potential at 
Tmf (i- e. at the bottom of the random band) is finite, the Tmf corresponding to the BEC 
would be zero |]l2l, due to the disordering effect of a large number of zero-energy (Goldstone- 
like) modes. Essentially, Tmf in this situation could vanish for the same reason as in the 
non-interacting gas of bosons in two dimensions. Motivated by this scenario I study first an 
admittedly artificial, but an instructive case of disordered lattice bosons with the infinite- 
range hopping, for which the mean- field theory is exact. In the infinite-range model the gap 
is uniform in space, and with decrease of particle density the transition temperature and 
the chemical potential indeed exhibit a clear BCS-BEC crossover if disorder is weak enough. 
Interestingly, Tmf nevertheless always remains finite, and the answer to the question posed 
above is negative. In the limit of low-density and week disorder the chemical potential drops 
below the band of random energies, and right at Tmf approaches the energy of the single, 
extended "bound" state, which then begins to become macroscopically occupied. The gap 
between this "bound" state and the rest of the band is finite at any disorder only for a 
uniform distribution of random energies; it collapses at strong disorder for any distribution 
that vanishes at the bottom, in which case the crossover to BEC of dilute bosons ceases to 
exist. Instead, at strong disorder one has the "localized superconductor" ground state [Q, 
where the superfluid state is phase-coherently formed out of macroscopic number of localized 
states. 

Retaining only the nearest-neighbor hopping in the effective bosonic theory leads to a 
more realistic and a more complicated problem, where the superconducting gap becomes 
strongly non-uniform in space. The non-linear mean-field equations that determine the 
values of the gap at different points in space can then be solved only numerically, however, 
some general statements on the nature of the superconducting solution can be made without 
explicitly finding the solution. Most importantly, the solution for the gap is either zero 
everywhere or nowhere in space. I discuss the mechanism behind the formation of the 
superconducting solution below Tmf, and the special role played by the sites with energies 
close to the chemical potential, which act as "seeds" of the long-range order 0]. In the limit 
of strong disorder these sites become sparse, and the gap exponentially quickly diminishes 
away from them, so that the system starts to resemble a granular superconductor with weak 
links connecting different granules. The transition temperature Tmf is determined by the 
seed which is in the best local environment to first trigger the long-range order. Within the 
mean-field picture Tmf is again always finite, and at T = amplitude of the superconducting 
gap is everywhere positive, although with diminishing average value as disorder is increased. 
I also comment on the additional separation of energy scales which occurs at strong disorder, 
where the superconducting gap and the minimum single-fermion energy cease to be the same. 



The paper is organized as follows. In the next section the effective low-energy bosonic 
Hamiltonian is derived starting from the disordered "negative-U" Hubbard model for elec- 
trons in the strong disorder limit. In Section 3 the BCS-BEC crossover at weak disorder and 
the "localized superconductor" at strong disorder in the model with infinite range hopping 
are discussed. The Tmf in the case of nearest-neighbor hopping and strong disorder is calcu- 
lated in the Section 4. Summary and connections to other works is given in the concluding 
section. 

2 The low-energy bosonic theory 

We will be interested in the system of lattice fermions in a random chemical potential, 
interacting via on-site attraction, defined with the standard negative-U Hubbard model: 

H = -t Yl 4,aCj,a+ Yl ViHi^^-UY^i^+^i,- (1) 

(j,,j)o-=+,— i,u=+,- i 

Here V^ is a random variable uniformly distributed in the interval [—W/4, W/4], and U > 0. 
When W/t is small, the eigenstates of the single-particle part of the Hamiltonian, given by 
the first two terms in Eq. 1, are close to Bloch waves, and the system has the usual BCS 
instability at small attraction U. In what follows I consider the opposite limit of a very strong 
disorder, when W/t >> 1. In the eigenbasis of its single-particle part the Hamiltonian (1) 
for a fixed realization of the random potential becomes: 

H = Y'^ana,a-Y< '^''^P\'^^'^> '^a,+na-- Y < «) TI^^I/?, ^ > cl,+C^,+4 -C^ _. (2) 

In the strong disorder limit the exact eigenstates of the random single-particle part of the 
Hamiltonian approach the site-localized wave-functions, \a >— ;> \i >, so < a,a\U\a,a >~ U, 
and the last term in the Eq. (2) which involves the overlap between different strongly 
localized states may be treated as a small perturbation. In the extreme site-localized limit 
when t/W = the last term in Eq. (2) vanishes, and the ground state is simply 

1^0 >= n 4-4+10 >. (3) 

i,{2Vi-U<0) 

Crucial observation is that the low-energy excitations above this ground state are purely 
bosonic: to create a pair of electrons at the next available energy level does not cost any 
energy, while to create a single electron there costs a finite amount of U/2 [Q. Despite 
the single-particle gap, the ground state does not have any off-diagonal long-range order 
characteristic of a superfiuid, since all the states are either fully occupied with a pair of 
electrons or empty, and it, in fact, represents a locahzed insulator. 

At strong disorder the last term in the Hamiltonian (2) is small, but finite. There are 
three different combinations of indices in this term that then give a dominant contribution: 
1) the Hartree term, a = P ^ •y = 6, 2) the pair-breaking term, a = S ^ (3 = •y, and 3) the 
pair "hopping" term, a = y ^ P = 6. 



For a zero-range interaction all three matrix elements are equal, and in the limit of strong 
disorder exponentially small compared to U. In contrast to the Hartree and the hopping 
terms 1) and 3), the term 2) annihilates any state that consists only of pairs of electrons, 
and therefore may be dropped in the low-energy sector of the theory, at energies much lower 
than U. One may then introduce the composite bosonic operators 

^o Cq.^-I_Cq.^_, Oq, Cq, _Cq, _|_ y^) 

and restrict the available Hilbert space by introducing the hard-core constraint as (6^)^ = 0, 
to account for the Pauli principle. The effective low-energy Hamiltonian in the strong- 
disorder limit is therefore: 

H = Y^Ejj)^ - Y. JaAbih + HKblbp), (5) 

with Ea = 2ea — U, and Ja,i3 =< a, a|t/|/3, /? >=< a, P\U\(y, P >■ In the strong-disorder limit 
we can further simplify the effective Hamiltonian (5) as follows. First, the random lattice 
formed out of maxima of sharply peaked localized states |a > is not very different from the 
regular lattice we started with, and we will assume that they are identical. The difference can 
be absorbed into the overlap integrals J^^^. Second, the random matrix elements Ja,i3 will be 
assumed to be all the same, and equal to J if the "sites" a and f3 are "nearest neighbors", 
and to zero if they are not. As will be shown later, this does not affect qualitatively the main 
features of the mean-field theory. The only random variables left in the effective theory then 
are the site-energies E^, and the effective Hamiltonian describes the system of hard-core 
bosons hoping on a regular lattice in presence of random chemical potential, i. e. the system 
of " dirty bosons" . 

The reader may have noticed that the transformation of the Hamiltonian from Eq. 1 
into Eq. 2 is completely general, and in principle may be performed even in the limit of 
weak disorder. In that case however, all the interaction terms in the Eq. 2 are of the same 
order, since the exact eigenstates are extended (or weakly localized in two dimensions) and 
the overlaps are numerous and large. So at weak disorder there is no clear separation of 
the Hamiltonian into the main part and the small perturbation, and the utility of the exact 
eigenstate representation is somewhat diminished. Proceeding nevertheless with the reduced 
bosonic Hamiltonian and assuming equal overlaps between the each pair of exact eigenstates 
leads to the Anderson theorem ||l|]. 

3 Infinite-range hopping 

It is often useful to consider first the modification of the theory for which the mean-field 
approximation should be exact. For that purpose, assume that the hopping integral Ja,i3 
in (5) is the same for all pairs {a, (3) on the lattice. As discussed above, this assumption 
would be justified if disorder in the initial Hamiltonian (1) were weak and the system is 
three dimensional, so that all the eigenstates \a > were extended in space: the overlap 



integral Ja^p would then indeed be nearly independent of the states a and (3. In the strongly 
disordered limit under consideration this is just a mathematical device. The effective low- 
energy Hamiltonian becomes 

H = Y.EXK-^Y.€h,, (6) 

where the second sum now is performed over all states a and /?, M is the total number of 
states (or lattice sites), and I dropped the Hartree term, which in the infinite-range model 
contributes only to an overall shift of the chemical potential. The hopping term is now of the 
form O'^O, with the operator O = J2a ba- Since (O) is an extensive quantity of order M, the 
standard mean-field decoupling O'^O -^ 0'^{0) + 0{0'^) — |(0)^| is exact, since the neglected 
fluctuation part is small (~ 1/Af) in the thermodynamic limit. The grand-canonical free 
energy for the infinite-range model is therefore given by [|l^ 



F = -^ -TY. InTrexp --[{E^ - Mb^ + mbl + m*6«], (7) 

where m = {J2aba)/J, and the chemical potential is determined by the total number of 
particles as A^ = —dF/dfi. 

The hard-core repulsion between bosons restricts the Hilbert space for the single site a 
to only two states, and the above trace can be easily evaluated. Up to constant terms, the 
disorder-averaged free energy per site is 



F Imp u r J(S-/x)2 + 4|m|2 

' ' ^ T P{E)dE In cosh y- 4- ^, (8) 



M J 2 J ' ' 2T 

where P{E) is the normalized distribution of random energies Ea- We will first take it 
to be constant for energies in the interval [—W/2, W/2] and zero otherwise; this should be 
adequate in the strong-disorder limit W/t » 1 of the electronic Hamiltonian (1). The order 
parameter m ^ and the chemical potential n are determined from the equations: 

1 = J P{E)dE^^^^^^, (9) 

•^ y'(E - /i)2 + 4|m|2 

and 



2T 

M 2 2J''~''"~ .>/(E-^)2 + 4|m|2 



^ ' ' ^ PiE)dE^^^P^^^^. (10) 



The transition temperature is determined by demanding that m = in the last two 

equations, so that 

J rW/2 tanh(,^) ^ ^ 

1 = — / dE \2Imf1 11 

wJ-w/2 E-fi ^ ' 



and 

Nil fW/2 E - a 

^. = --^77 dEtanh(=— -^). (12) 

M 2 2W J-w/2 2Tmf 

These two equations in general need to be solved numerically; there are two simple limits 
however, which allow a transparent analytic solution. First, consider the density of bosons 
to be N/M ^ 1/2. As the solution for Tmf in this case is finite, it immediately follows from 
the second equation that ^{Tmf) ~ 0, i. e. near the center of the band of random energies. 
More precisely, in the strong-disorder limit when J/W « 1, niTMp) ~ W{N/M — 1/2) 
near half-filling. The Eq. (11) then implies that 




N W 

MF ^const.PF^/ — (1-— )exp(-— ), (13) 

which has the standard BCS form. In the strong disorder limit under consideration J « W, 
and consequently T^jp is much smaller than in the clean limit. As the particle density 
decreases from half-filling ^{Tmf) becomes negative, and Tmf decreases from its maximum 
value. At small density the simple linear relation between the chemical potential and the 
density of bosons breaks down, as the system crosses over to a dilute limit of almost non- 
interacting bosons. At exactly A^ = we see from the Eq. (12) that Tmf = and ii{Tmf) < 
-W/2. It follows from Eq. (11) that as A^ ^ 

W W 

MrMF)--yC0th(— ) (14) 

from above. The critical temperature when A^ — > may be obtained approximately: 

iyexp(-iy/J) 



■MF 



■ln(A^/M) 



(15) 



it logarithmically slowly goes to zero with decreasing density. 

In the dilute limit N/M — > 0, the above results for the transition temperature and the 
chemical potential become identical to those at the BEC of non-interacting bosons in the 
same random potential. If one would neglect the hard-core constraint on the occupation 
number of bosons, the Hamiltonian (6) would describe non-interacting particles, and its 
disorder-averaged spectrum could be calculated. The average density of eigenvalues in the 
thermodynamic limit is []T1| 



p(A) = P(A) + ^5(A + Ao), (16) 

where -P(A) is exactly the same distribution of random energies as in Eq. (8), and — Aq < 
—W/2 is the energy of the single extended bound state given by fl^ 



^^V^^aTv- ("' 



A + Ao 



This is nothing but the Eq. (11) with fi replaced by — Aq, in the hmit when Tmf ~^ 0. 
Thus, — Aq = ii{Tmf)^ as in Eq. (14). The temperature of BEC into this bound state for the 
non-interacting bosons would be given by 

" 'pw .it , . (18) 



M J 'exp(^)-l 

which differs from the Eq. (12) only in sign in front of unity in the denominator. This 
becomes negligible in the limit of small number of bosons when Tbec is small, and as a 
result, the transition temperature calculated from the Eq. (11) approaches the temperature 
of BEC for the non-interacting system. 

The reason for finiteness of Tmf even in the dilute limit in the infinite-range model lies 
in the fact that the energy of the bound state is always below the rest of the band, at least in 
the case of uniform distribution of random energies studied above. Pazmandi et al. [|1^ have 



also shown that in this case the bound state is spatially extended, and is given by a linear 
combination of macroscopically large number of localized states \a >. This is the reason 
why it is able to accommodate a macroscopic number of hard-core bosons. This breaks 
down for distributions that tend to zero at the lower edge; for example, for the triangular 
distribution P{E) = 2{E + W/2)/W^ for -W/2 < E < 0, and P{E) = 2{W/2 - E)/W^ 
ioT < E < W/2, it follows from Eq. (11) that as Tmf ~^ 0, /^(Tmf) < -W/2 only 
for W < (2 In 2) J (see Figure 2). For disorder stronger that this critical value chemical 
potential at the transition temperature lies in the band of localized state, and there is no 
condensation into any particular state. This reflects the localized nature of the states for 
this distribution, none of which, unlike the bound state in the weak-disorder limit, can 
accommodate a macroscopic number of hard-core bosons. We see that the case of uniform 
distribution (or more precisely, of distribution which remains finite at the lower edge) is 
indeed quite special; in general, the crossover to BEC in the dilute limit happens only at 
weak disorder. At stronger disorder the transition at Tmf does not correspond to BEC into 
any particular state, but to a macroscopic number of localized states forming a condensate 
with a well defined phase. This is precisely the state Ma and Lee described as a "localized 
superconductor" ||^. 

4 Nearest-neighbor hopping 

I now turn to a more realistic case of nearest-neighbor hoping of bosons and assume a fixed 
chemical potential at /x = 0, i. e. half-filled band. The mean-field free-energy in a fixed 
realization of disorder is now 



MF 



_g2 _j_ M^fji |2 

J2 ^*a Ja^p^P -Tj2^n cosh ^ " — —, (19) 



2T 



where the site-dependent order-parameter rria is determined as 



rrir 



mp 



Ei + 4\mp 



2 



JV , " tanh^^ — . (20) 



/3 



^| + 4|m^|2 2T 



The Hartree term is again neglected, but this time purely for convenience; its presence would 
merely make the analysis more complicated without affecting the main points. The sum in 
the last equation is performed only over the nearest neighbors of the site a. 

One faces a set of M coupled non-linear equations, and at first it is not obvious that 
a non-trivial (ttIq, ^ 0) superconducting solution exists at all. One can notice from the Eq. 
20 however, that if it does, all ma can be chosen to be real. It is useful to consider a toy 
two-site problem. At T = the solution of the gap equations is 



"^1(2) = J^ 



1 - {E,E,/J^f 



The real solution exists only if E1E2/ J"^ < 1, i. e. if the random energy of at least one 
of the sites is sufficiently close to the chemical potential. In a lattice there obviously will 
be many nearest neighbors where this is not satisfied; nevertheless, the superconducting 
solution always exists due to the lattice sites with energies close to /x, which then spread 
the long-range order throughout the whole lattice 0]. To see how this happens consider a 
lattice where all sites have energies in vicinity of some large energy E = Vt » J , so that 
the only possible solution of the gap equations is the trivial {iria = 0) solution. Then change 
the energy of a single site at a = into Eq = 0. The solution at T = in the strong disorder 
limit to the lowest order in J/Q is: 



mo = 
at the site a = 0, and 



•"» = IF + '''??>• '22) 



"■> = •"■Ti''"' + '^'•''^*"* '23) 

at all other sites, where integer n measures the shortest lattice distance of the site a from 
the site (n = 1 is the nearest neighbor), z is the lattice coordination number. Thus a single 
site with the energy at the chemical potential "polarizes" immediately the whole lattice at 
T = 0, and the gap opens everywhere. The same mechanism is responsible for the existence 
of the non-trivial solution in a truly disordered lattice with a distribution of energies. In the 
thermodynamic limit, for any distribution that is finite arbitrarily close to the Fermi level 
(i. e., does not have a true gap there) there always will be sites with energy arbitrarily close 
to fi. At large disorder, the solution of the gap equation will be qualitatively the same as 
in the above model of a two-component alloy, where the majority of sites has a large energy 
~ W, and there are rare "seeds" with zero energy in between. Allowing the values of J to 
vary in space can not affect the conclusion either; since hoping integrals are all necessarily 



positive, even if they are allowed to be random does not change the fact that the non-trivial 
solution of the gap-equation is positive everywhere. The ground state will therefore always 
be a superfluid within the mean-field theory. 

Having this in mind it becomes a simple matter to obtain the critical temperature at 
which the non-trivial solution first becomes possible. Consider such a "seed" of long-range 
order at a = and its nearest neighbors; their energies will typically be of the order of W/2. 
Since Tmf ~ J, in the limit J/W << 1 the equations become: 

^o-JE/^' (24) 

where I assumed that both Tmf and rrifj << l-E/j], and 

771(3^ J tanh—^. (25) 

Thus the gap at the site a = satisfies the self-consistent equation 

^o^j'E^tanh^, (26) 

which allows a non-trivial solution below the critical temperature 

1 






In a given realization of the random potential Tmf will be determined by the "seed" for 
which the above expression is maximal. Typical value of the mean-field critical temperature 
is thus 

Tmf ^ ^, (28) 

roughly equal to the value of the gap at the site a = at T = 0, and quite unlike the 
BCS-like expression (13) in the infinite-range model. Most importantly, Tmf is again always 
finite, and vanishes only in the limit of infinite disorder. This is in accordance with the result 
that the ground state within the mean- field theory is always superconducting [Q. 

5 Discussion 

We argued that in the negative-U strongly disordered Hubbard model the effective theory at 
energies much below U is completely bosonic, and then proceeded to obtain the mean-field 
transition temperature at which the bosonic electron pairs would enter the superfluid state. 
True Tc is always below our mean- field result, and in particular, the quantum fluctuations 
of the phase neglected in the mean- field theory should drive it to zero at a finite disorder. 
The obtained Tmf should be understood as a crossover temperature where the superfluid 

9 



order parameter begins to develop an appreciable amplitude. Our main conclusion is that 
the Tmf is always finite, and vanishes only in the limit of infinite disorder. 

In the pure system and at weak coupling, Tmf, or the zero-temperature superconducting 
gap, also measures the binding energy of an electron pair, i. e. the gap for single-electron 
excitations. In the strongly disordered system studied here, localization of the single particle 
wave-functions makes the pairs tightly bound, and it is f/ >> Tmf that determines their 
binding energy. In the mean- field ground state at T = 0, the superconducting gap and 
the minimum energy of single-fermion excitations, while close at weak disorder, become 
completely different at stronger disorder (see Figure 3). Initially both are suppressed at 
weak disorder essentially because of the decrease in density of states at the Fermi level; 
however, while the spatial average of the superconducting gap continues to decrease with 
increasing disorder, single-fermion gap after initial suppression must begin to increase, since 
at very strong disorder it approaches U. The turning point marks the entrance into the 
strong-disorder limit. This divorce of the two energy scales due to disorder is clearly seen in 



numerical calculation of Ghosal et al. |T5|, where the spectral gap of Bogoliubov - de Gennes 
excitations (which correspond to breaking of pairs) starts increasing at large disorder. This 
is quite similar to what happens in the clean model with increase of attraction between 



fermions |TT| . In sum, there are three characteristic energy scales in the disordered problem: 
the pair binding energy ~ U, the zero-temperature superconducting gap ~ Tmf, and the 
true Tc where the phase coherence sets in, and which lies beyond the scope of the mean-field 
considerations. 

In the simplified effective theory of the hard-core disordered bosons with infinite range 
hopping and with a bound distribution of random energies that vanishes at the bottom, 
there is a crossover to the BEG in the dilute limit if disorder is weak enough. At stronger 
disorder one is in the BGS "localized superconductor" limit [Q, where a large number of 
localized states forms a condensate with a well defined macroscopic phase. In both cases 
the mean-field transition temperature remains finite at any disorder. There is again a close 
analogy between these results and the BGS-BEG crossover in the clean system of attracting 



fermions when the strength of attraction is varied [|Tl|]. In that case the role of hopping 
J is played by the attractive interaction, and that of disorder W by the band width t: in 
three dimensions density of states at the bottom of the band vanishes, and the bound state 
exist only if attractive interaction is large enough. As one increases the attraction at a fixed 
density of particles the system crosses from the BGS to the BEG limit, where the bound 
state below the band becomes macroscopically occupied. The case of uniform distribution of 
random energies in our problem would in this sense be similar to the case of two dimensions, 
where the bound state exists already for an infinitesimal attractive interaction. 

It is interesting that although at first sight the nature of the BEG and that of the 
localized superconductor is quite different, the crossover between them in the infinite range 
model is completely smooth, without a phase transition. It has been argued that in the 
strong disorder limit, at T = 0, if one enters the superfiuid phase by varying chemical 
potential, when /i is precisely at the lower edge of the band system is in the Bose-glass phase 



Q . According to this scenario one could have a direct Mott insulator - superfiuid (of BEG 

10 



type) transition at weak disorder, and Mott insulator - Bose glass - superfluid (localized 
superconductor) sequence of transitions at strong disorder, depending whether there is, or 
there is not, the extended bound state in the spectrum. From our point of view the alleged 
Bose glass phase corresponds to the zero density of bosons, or more pictorially, to just a 
few bosons localized in the random potential. It is therefore not a macroscopic phase. At 
any finite density the mean-field theory yields a superfluid ground state of either BEC or 
localized superconductor type, with a smooth crossover between them. 

The mean-field transition temperature in the model with only nearest-neighbor hopping 
is also always finite, and in the strong-disorder limit may be obtained by considering the 
opening of the superconducting gap at a particular zero-energy site. This immediately 
triggers the long-range order in the entire lattice, and the gap opens everywhere below Tmf- 
In the strong disorder limit superconducting gap vanishes exponentially fast away from these 
special zero-energy sites. As disorder is increased, both Tmf and the average value of the gap 
decrease, but always remain finite. The low energy excitations are located at the zero-energy 
sites, where the gap is ~ J'^/W at T = 0, while at farther sites the value of the gap is very 
close to zero. This has been observed in the numerical solution of Bogoliubov-de Gennes 



equations by Ghosal et al. [|T5[, who find spatial correlations between the "superconducting 
islands" where the gap is substantial and the sites with small random energies. At large 
disorder the system becomes very much like a granular superconductor. This feature can be 
readily understood in the strong-disorder limit. 

Evidently, the description of the superfluid-insulator quantum phase transition in the 
simple attractive model studied here requires an approach beyond the mean-field theory [Q, 
1^, []^. The mean-field theory does indicate however, that the relevant degrees of freedom for 
the long wave-length physics near such a transition are quantum phase fluctuations, which 
become enhanced at strong disorder due to a small typical amplitude of the order parameter. 
In ref. 7 it was shown that a disordered Bose system will turn into an insulator at T = 
if the average phase stiffness at a microscopic scale is small enough. Small randomness in 
the stiffness, however, turns out to be irrelevant at the superfluid-Bose glass transition, at 
least within the e = d — 1 expansion [^ . Thus the model of interacting bosons in random 
potential should be the appropriate starting point for addressing the universal properties of 
the superconductor- insulator transition, at least in the strong disorder limit. We may obtain 
a crude estimate of the shape of the superfiuid-Bose glass phase boundary at strong disorder 
as follows: dirty-boson Hamiltonian in Eq. 5 will have an insulating Bose-glass state at some 
critical disorder W = xJ, where a; is a non-universal number of order of unity. Since the 
random eigenstates are strongly localized, away from its maxima we may assume they go 
to zero exponentially, over a characteristic length (in units of lattice spacing) oc t/W. The 
overlap integral then J oc exp—iW/t)"^, so that at strong disorder the superfiuid-Bose glass 
phase boundary is at: 

-(xexpc(— )2 (29) 

where c is yet another non-universal constant. In principle there could be a prefactor pro- 
portional to some power of VT/t, but the main point is that at strong disorder one expects 
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the transition to occur at a very large U. Thus the phase boundary in the U — W plane 
when observed over a limited range of large U should resemble a constant, in qualitative 
agreement with the recent Monte Carlo calculation [p!7[] . 

The crucial question is if the purely bosonic nature of the superconductor-insulator 
transition persists at arbitrarily small U. It has been claimed recently [|I^ that at small 
attraction the exit from the superconducting state coincides with the appearance of the 
zero-energy single-fermion excitations, i. e. with breaking of electron pairs. I suspect that 
this feature may be a finite-T artifact of the numerical simulations, in which temperature 
enters as a finite length of the system in the imaginary time direction. In fact, when t = in 
the Hamiltonian (1) the single-particle density of states can be calculated exactly |]16|, and 



while at T = there is always a gap equal to U, ii U < W/2 density of states at zero energy 
becomes finite at any finite temperature. The finite temperature thus can mimic the closing 
of the single-fermion gap if the attractive interaction is weak enough, even though at zero 
temperature the gap is always present. 

Another important question is validity of the bosonic theory if repulsion between elec- 
trons is included. Increasing disorder localizes the wave-functions and thus decreases the 
"size" of the electron pairs; this effectively increases the repulsion and may lead to pair- 
breaking. A precursor of this effect is the reduced screening of the Coulomb interaction due 
to disorder, which is believed to be an important source of reduction of order parameter 
amplitude or Tc in real systems P]. It is conceivable that if the repulsion is large enough pair 
breaking becomes the mechanism of the transition, so that the transition out of the super- 
fluid is into an electronic glass. Another possibility is that the Bose-glass state will always 
set in due to quantum fluctuations when the average amplitude of the superconducting order 
parameter becomes sufficiently small, so that the loss of phase coherence always occurs before 
pairs dissociate. In this scenario, at strong repulsion Bose-glass would exist only in a narrow 
region between the superconductor and the electronic glass, where the electron pairs become 
finally broken. Studies of the interplay of the two effects are obviously required before our 
understanding of the superconductor-insulator transition may be considered complete. 

6 Acknowledgment 

This work has been supported by NSERC and the Faculty of Science research grant D-5765 
at Dalhousie University. Part of it was also supported by Killam foundation while the author 
was at the University of British Columbia. 



12 



FIGURE CAPTIONS: 



Figure 1: The schematic temperature-disorder phase diagram of the negative-[/ dis- 
ordered Hubbard model. While the Tmf (dashed line) initially depends on disorder only 
weakly (the regime of Anderson theorem), for stronger disorder it becomes suppressed (Ma 
and Lee regime). Tmf always remains finite at finite disorder. True critical temperature T^ 
(full line) vanishes entirely due to quantum phase fluctuations. The thick line at T = and 
W > Wc represents the Bose glass insulator. 

Figure 2: Example of the density of states in the infinite-range effective theory below 
(a) and above (b) critical disorder |1^. Below, there is a single eigenstate (weight 1/M) 



detached from the band, which is spatially extended, and available for Bose condensation 
at low densities. The rest of the band is localized. Above the critical disorder the extended 
state merges with the continuum and localizes. Nevertheless, the ground state of the system 
is always a superconductor. 

Figure 3: Qualitative behavior of the average superconducting gap at T = (A ~ Tmf) 
and the single-fermion gap (dashed line) with disorder. While the two are close at weak 
disorder (and not too strong U), at large disorder they differ by U. 
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